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ABSTRACT 

We combine information on cosmological parameters from cluster abundances, CMB primor- 
dial anisotropics and the IRAS 1 .2 Jy galaxy redshift survey. We take as free parameters the 
present values of the total matter density of the universe, ilm, the Hubble parameter, h, the lin- 
ear theory rms fluctuations in the matter density within 8/;^' Mpc spheres, Og, and the IRAS 
biasing factor, ^jiras- We assume that the universe is spatially flat, with a cosmological con- 
stant, and that structure formed from adiabatic initial fluctuations with a Harrison-Zel'dovich 
power spectrum (i.e. the primordial spectral index n~\). The nucleosynthesis value for the 
baryonic matter density £2b = 0.019//;^ is adopted. We use the full three- and four- dimen- 
sional likelihood functions for each data set and marginalise these to two- and one- dimen- 
sional distributions in a Bayesian way, integrating over the other parameters. It is shown that 
the three data sets are in excellent agreement, with a best fit point of = 1 — = 0.36, 
h = 0.54, ag = 0.74, and /jjras = 1.08. This point is within one sigma of the minimum for 
each data set alone. Pairs of these data sets have their degeneracies in sufficiently differ- 
ent directions that using only two data sets at a time is sufficient to place good constraints 
on the cosmological parameters. We show that the results from each of the three possible 
pairings of the data are also in good agreement. Finally, we combine all three data sets to 
obtain marginalised 68 per cent confidence intervals of 0.30 < £lm < 0.43, 0.48 < h < 0.59, 
0.69 < Og < 0.79, and 1.01 < ^^iras < 116. For the best fit parameters the CMB quadrupole 
is 2rms-ps 18.0/jK, the shape parameter of the mass power-spectrum is F = 0.15, the baryon 
density is D.\, = 0.066 and the age of the universe is 16.7 Gyr. 

Key words: large-scale structure of Universe - cosmic microwave background - galaxies: 
clusters 



1 INTRODUCTION 

Combining and comparing information from different types of ob- 
servations is a powerful and important tool in cosmology. In this 
letter we investigate the agreement between the conclusions drawn 
from studies of CMB data, a galaxy redshift survey and the varia- 
tion of the cluster abundance with redshift. 

The number density of rich galaxy clusters at low redshift de- 
pends strongly on both and Cg, with a weak dependence on 
the shape of the power spectrum of fluctuations (Lilje 1992; Bah- 
call & Cen 1993; Hanami 1993; White, Efstathiou & Frenk 1993). 
Measuring the evolution of the number density with redshift breaks 
the degeneracy between these two main parameters (eg. Oukbir & 
Blanchard 1992; Viana & Liddle 1996; Eke, Cole & Frenk 1996). 
The comparison of the X-ray cluster temperature functions (the 
number of clusters per unit volume as a function of the temper- 



ature of their X-ray emitting gas) determined at low (Edge et al. 
1990; Henry & Amaud 1991) and high (Henry 1997) redshift is a 
good way to implement this test. Henry (1997) first performed this 
operation, and Eke et al. (1998) and Viana & Liddle (1999) have 
considered additional systematic uncertainties. 

The power spectrum of the primary anisotropies in the Cos- 
mic Microwave Background (CMB) depends on many cosmolog- 
ical parameters. The constraints that can be inferred from current 
CMB data arise mainly from the amplitude measured for the fluc- 
tuations on large angular scales and also the position and height 
of a peak in the power spectrum at smaller scales. Discussions of 
the implied cosmological parameter information, given a range of 
assumptions, have been most recently presented by Hancock et al. 
(1998), Bond & Jaffe (1998), Lineweaver (1998), Tegmark (1998), 
Bartlett et al. (1998) and Efstathiou et al. (1998). 
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While the CMB and cluster abundances depend on the fluc- 
tuations in mass, galaxy redshift surveys tell us about the distri- 
bution of luminous matter. It is commonly (and probably naively) 
assumed that the distribution of luminous galaxies of a certain type 
is related to the underlying matter distribution by the 'biasing fac- 
tor', b. We define this as the ratio of the rms fluctuations in the 
galaxy number density to those in the underlying mass. Thus, the 
amplitude of the galaxy power spectrum is larger than that of 
the mass. Under the assumption that b is scale-independent, the 
shape of the galaxy power spectrum is the same as that of the mat- 
ter power spectrum. Thus the galaxy power spectrum can be used 
to constrain the matter power spectrum shape parameter, T. As the 
cluster abundances only weakly constrain F, the addition of galaxy 
redshift surveys to the joint analysis helps to fix the shape of the 
power spectrum. In addition, the redshift space distortion of the 
galaxy distribution due to peculiar velocities yields a measure of the 
combination, (3 = flUj^/Z?. The joint analysis with CMB and cluster 
abundances allows us to remove the degeneracy in this combina- 
tion. Constraints on cosmological parameters from redshift surveys 
have been obtained by Fisher, Scharf & Lahav (1994), Tadros et al. 
(1999), Sutheriand et al. (1999) and Branchini et al. (1999). 

Many authors have now compared various combinations of 
data sets at a range of levels of detail (Bond & Jaffe 1998; Gawiser 
& Silk 1998; Lineweaver 1998; White 1998; Garnavich et al. 1998; 
Webster et al. 1998; Eisenstein, Hu & Tegmark 1998; Efstathiou 
et al. 1998 and Roos & Harun-or-Rashid 1999). We choose to per- 
form a joint analysis of CMB, galaxy redshift survey and cluster 
abundance data sets, since they provide complementary constraints 
on cosmological parameters. As we demonstrate below, the param- 
eter degeneracies arising from these three data sets are relatively 
complicated and merit the full likelihood approach adopted in this 
paper. 

In Section 2 we introduce the data and outline the methods 
used to produce likelihoods for each separate set of observations. 
We describe our method of combination and some statistical issues 
in Section 3 and present our results in Section 4, discussing them 
in more detail in Section 5. 



2 INDIVIDUAL DATA SETS 
2.1 Cluster Abundances 

A range of different values for cosmological parameters have been 
inferred from the evolution of the number density of clusters with 
redshift (e.g. Frenk et al. 1990; Carlberg et al. 1997; Bahcall, Fan & 
Cen 1997; Blanchard & Bartlett 1998; Reichart et al. 1999), and in 
some cases the quoted uncertainties are sufficiently small to make 
the results incompatible. This implies that systematic uncertainties 
are important (see Eke et al. 1998 for a discussion). For a joint 
likelihood analysis, the resulting confidence limits on the inferred 
cosmological parameters are only reliable if the systematic errors 
present in the individual analyses are less important than the statis- 
tical ones. For this reason, the cluster abundance is addressed using 
the X-ray temperature function, where this is least unlikely to be 
the case. 

We use two X-ray flux-limited samples. The first is comprised 
of the 25 clusters with average redshift 0.05 compiled by Henry & 
Arnaud (1991). We have incorporated the latest temperature mea- 
surements from Ginga and ASCA for these clusters (Henry 1999). 
The second sample is comprised of 14 Einstein Extended Medium 
Sensitivity Survey (EMSS) clusters with average redshift 0.38 com- 
piled by Henry (1999). The temperatures of all 14 clusters have 



been determined with ASCA data using the most recent calibrations 
available (Henry 1999). In order to calculate likelihoods the Press- 
Schechter (1974) expression for the abundance of rich clusters as 
a function of mass and redshift, coupled with a mass-temperature 
conversion (White et al. 1993; Eke et al. 1998) is employed to pre- 
dict the distribution of cluster redshifts and temperatures for each 
model. These are compared with the observed distributions as dis- 
cussed by Henry (1997) and Eke et al. (1998), using the 'default' 
assumptions detailed by Eke et al. (1998). 

As discussed in Section 1, the main constraint from cluster 
abundances is in the Og, flm plane, and this is plotted in Fig. |l| (a). 

2.2 CMB 

Our approach follows that of Hancock et al. (1998) and Webster et 
al. (1998). We use the same data as Webster et al., with the addition 
of the recent QMAP points (de Oliveira-Costa et al. 1998). This 
same data set is used and plotted in Efstathiou et al. (1998). 

The flat band-power method is employed (as described in e.g. 
Hancock et al. 1998), although Bartlett et al. (1999) and Bond, 
Jaffe & Knox (1998) find that this approximation widens the er- 
ror ranges on the cosmological parameters. The conclusions of this 
letter would not be significantly affected by such a change. 

We assume that there are negligible tensor contributions to the 
CMB power spectrum, as predicted by most inflation models. We 
also assume neglible re-ionisation and a Harrison-Zel'dovich pri- 
mordial power spectrum (ie. the primordial spectral index n = \) 
as favoured by inflation and the CMB data. We note that chang- 
ing these assumptions would substantially alter our results. The re- 
maining family of Cold Dark Matter models are analysed using the 
Boltzmann code of Seljak & Zaldarriaga (1996). 

The CMB COBE data point constrains the large scale temper- 
ature fluctuations well. This converts to a strong constraint on 
for each h and Q.^. Estimates of the location and height of the first 
acoustic peak are also provided by the CMB data, and these con- 
strain a degenerate combination of h and Sl^ (Fig. |l| (b)). Putting 
these together also places a weak constraint on Og. 

2.3 IRAS 1.2 Jy survey 

The data and method used here are exactly as in Fisher et al. (1994) 
and Webster et al. (1998). We use a sample of 5313 IRAS galaxies 
flux-limited to 1.2 Jy at 60 /jm, selected from the IRAS data base 
(Strauss et al. 1990, Fisher 1992), and calculate redshift weighted 
spherical harmonics which may be compared with those predicted 
given a matter power spectrum. We use Gaussian windows for the 
redshift weighting, and take into account the redshift distortion us- 
ing an approximation that is correct only in linear theory. Therefore 
in this analysis we use only spherical harmonics up to / = 10 so that 
the scales probed are large, and the fluctuations can be assumed to 
be in the linear regime. 

It is assumed that the relationship between the galaxy distribu- 
tion and the matter distribution can be described by the single pa- 
rameter, the biasing factor for the IRAS galaxies, binA5 - In practice 
the relationship between the two power spectra may have many de- 
pendencies, including scale dependence. Note that if ^iras is actu- 
ally scale dependent then the information on the shape of the power 
spectrum given by the IRAS data would be wrong. However, if the 
assumption of scale independent biasing is correct then galaxy sur- 
vey data such as the IRAS 1.2 Jy survey place powerful constraints 
on the shape of the power spectrum in the scale range examined. 
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Figure 1. (a) The cluster abundance likelihood function in the {i2n,,0g}-plane, after marginalisation over h. (b) The CMB likelihood function in the {Q.^, 
/i}-plane, after marginalisation over Og. (c) The IRAS likelihood function in the {Q.^, ftj-plane, after marginalisation over Cg and ^iras- The contours denote 
the 68, 90, 95 and 99 per cent confidence regions. 



The redshift distortion also provides limits on Piras = f^m^/^lRAS 
and hence a constraint on Q.^ for each value of ^iras- However, 
this turns out to be a relatively weak constraint. Nevertheless, the 
range of values of Q.^ allowed by IRAS together with CMB or clus- 
ter abundances lies well within the range allowed by IRAS alone. 
Thus, for this paper, the most important constraint from the IRAS 
data is that on the shape of the matter power spectrum which, as 
will be seen in the next section, manifests itself in the fim, h plane, 
plotted in Fig. |^ (c). 



3 COMBINING THE DATA SETS 
3.1 Parameters 

In this letter the three cosmological parameters we vary are: 
the present day value of the reduced Hubble constant, h — 
Hq/{WO kms^'Mpc^' ), the matter density of the universe in units 
of the critical density, flm, and the matter power spectrum normal- 
isation parameter, Og, the rms linear fluctuation of matter density 
in 8/;^' Mpc spheres. The fourth parameter is ^iras. which should 
be regarded as a 'fudge factor' reflecting our ignorance on how to 
relate fluctuations in galaxy counts to fluctuations in mass. We as- 
sume that the universe is flat so that flm + Q./^ = 1 as suggested by 
recent analyses of CMB data which allow for flni and fl^ to vary 
independently (eg. Lineweaver 1998; Tegmark 1998; Efstathiou et 
al. 1998). 

To transform the information on the shape of the matter power 
spectrum provided by the IRAS and cluster abundance data into 
constraints in the parameter space we consider, we use the Bardeen 
et al. (1986) approximation formula for the matter power spectrum, 
parameterised by the shape parameter, P. We use the approximate 
expression provided by Sugiyama (1995), 



r = Q.mh exp — f2| 



fin 



(1) 



This, and the CMB power spectrum calculation, require a value for 
the baryon density which we take from the nucleosynthesis con- 
straint, ilb/!^ = 0.019 (Buries & Tytler 1998). Our results are rela- 
tively insensitive to this value. 

Using this conservative number of parameters allows a more 
thorough examination of the parameter space, and despite our large 
number of assumptions, we find remarkable agreement between the 
data sets. 



3.2 Statistical issues 

Comparing results from different data sets enables us to shed some 
light on how well our current picture of cosmology stands up. This 
kind of comparison could be performed by simply comparing the 
parameter values predicted by different data . While this is a start, 
it does not take into account the fact that most measurements are 
sensitive to more than one parameter. A full examination of the 
level of agreement between data sets should include a comparison 
of likelihood values in multi-dimensional space. In this letter we 
have considered the four free parameters flm, Og, h and ^iras and 
found that there is a region in the four-dimensional space in which 
the likelihood functions of each data set are high. This is important 
to check before simply multiplying likelihoods, because the final 
best fit point could be a bad fit to all the data sets used, even if the 
one-dimensional marginalised joint likelihood functions for each 
parameter look promising. The goodness of fit of each data set can 
be evaluated e.g. by x^. 

Combining results is also powerful because the parameter de- 
generacies from each data set may be complementary. Fig. ^ shows 
a three-dimensional view of the orthogonality of the CMB and clus- 
ter abundance confidence surfaces. The favoured region of parame- 
ter space is very much smaller than the regions of parameter space 
allowed by the single data sets alone, thus yielding much tighter 
constraints. This joint likelihood function is computed by calcu- 
lating likelihoods over a grid in parameter space, then multiplying 
them together for the different data sets. 

Note that multiplication of the likelihoods for each data set is 
strictly correct only in the case where the data sets are indepen- 
dent. While the CMB clearly probes a different part of the universe 
to either the IRAS galaxies or the cluster abundance clusters, it is 
not so clear whether there is some overlap in the IRAS and low 
redshift cluster data. In this letter we assume that these two data 
sets are independent, and justify it on the basis that (a) the redshift 
ranges do not overlap significantly, the low redshift clusters having 
a median redshift of 0.05, while the IRAS galaxies have a median 
redshift of 0.02 and (b) because IRAS galaxies are selected in the 
infra-red, they tend to be field galaxies and so are not likely to be in 
the clusters sampled in the cluster survey. This issue also affects the 
CMB analysis, since the flat band power analysis is not valid if 
the observations used overlap on the sky and in angular scale range 
(spherical harmonic I range). We have avoided this by using a set 
of CMB data points that do not overlap in this way, but the set of 
data points chosen could be argued to be somewhat subjective. 
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Figure 2. CMB and cluster abundance iso-probability surfaces plotted in 
three-dimensions, for Alog(Likelihood)=7.1 (3c for a 3d Gaussian proba- 
bility distribution). 

We obtain confidence limits on each cosmological param- 
eter by marginalising over the other parameters. If the parame- 
ters are well enough constrained then the integration limits used 
for the marginalisation will not affect the results. The Bayesian 
interpretation of the integration limits is that we have 'top hat' 
priors on each parameter. For the integration in the marginalisa- 
tion we use 0.2 < < 1.4, 0.3 <h< 0.9, 0.1 < fJm < 10 and 
0.7 < foiRAS < 2.0 (in 100, 100, 100, 50 steps respectively). The 
parameters are so poorly constrained by any single data set that 
these integration limits do affect some marginalised results. For ex- 
ample, our range of integration in Og affects the two-dimensional 
plot of h and Cl^ for the CMB data (Fig. |j| (b)), making the con- 
tours foreshorten at low Sl^. However, we choose our limits to be 
wide enough that they do not affect the results when using more 
than one data set at a time. Thus the one-dimensional marginalised 
distributions for pairs of data sets plotted in Fig. ^ (a) are inde- 
pendent of the priors. This integration process is computationally 
costly, and some authors utilise approximations at this point. For 
example Lineweaver (1998) and Tegmark (1998) 'project' the like- 
lihood function to one dimension, taking the maximum value of 
the likelihood found on varying the other parameters, rather than 
integrating over all the values of the likelihood for all the possible 
values of the other parameters. This approximation is correct when 
the likelihood function is a Gaussian, but not otherwise. Note that 
when we talk about the 68 per cent confidence surface or contour 
we mean the (iso-probability) boundary to the region containing 68 
per cent of the probability (given our top hat priors). 



4 RESULTS 

Since it is possible to choose reasonable integration ranges such 
that the results using any pair of data sets are independent of the in- 
tegration ranges, we may say that just a pair of data sets is sufficient 
to constrain properly these four parameters (or three parameters for 
CMB and cluster abundances). 



Parameter 


Best fit point 


68 per cent confidence limits 




0.36 


0.30 < Q.^ < 0.43 


h 


0.54 


0.48 < h < 0.59 




0.74 


0.69 < Og < 0.79 


i'lRAS 


1.08 


1.01 < foiRAs < 1.16 



Table 1. Parameter values at the joint optimum. The 68% confidence limits 
are shown, calculated for each parameter by marginalising the likelihood 
over the other variables. 

The results for pairs of data sets are shown in Fig. ^ It can 
be seen that these three data sets are in extremely good agreement, 
since the parameter values preferred by each pair are consistent. 
This is particularly apparent for Q.^. It is a result of the fact that 
all three 68 per cent confidence surfaces intersect in the same place 
in three and four dimensions. It was therefore not surprising to find 
that the best fit point given in Table 1 is within the 68 per cent 
confidence surface for each single data set. 

For the parameter values at the best fit point, the shape 
parameter for the mass power spectrum F = 0.15 ±0.02, the 
baryonic density flj, = 0.066 ± 0.013, the CMB quadrupole 
2rms-p.s = 18.0± I.O/jK, the IRAS distortion parameter Piras = 
fim''/^lRAS = 0.50 ± 0.05 and the age of the universe is 16.7 ± 1 .0 
Gyr. In each case, the errors correspond to the estimated 68 per cent 
marginalised confidence limits. 

As an indicator of how well the joint optimum fits each data 
set alone we calculate the value at the joint optimum for each of 
the CMB and IRAS data sets, using the method described in Web- 
ster et al. (1998). The CMB is 22.8 (we use 23 CMB data points) 
and the IRAS is 142 (for spherical harmonics up to / = 10 there 
are 120 degrees of freedom), suggesting that the predicted obser- 
vations agree well with the actual observations for the parameters 
at the joint optimum. Because of the unbinned nature of the cluster 
abundance analysis, it is more complicated to find a similar statis- 
tic. However, as discussed by Eke et al. (1998), all of the models 
considered predict a relatively even distribution of clusters over the 
range 0.3 < z < 0.4, but in actual fact there are many more clusters 
at the low end of this redshift range than predicted by any of the 
models considered. This suggests that the selection function for the 
clusters was poorly modelled, inevitably leading to a poor good- 
ness of fit to the cluster data, independent of the parameter values 
used in the model (given the assumptions of this paper). The num- 
ber of clusters predicted from the best fitting model is 19.5 at low 
redshift and 26 at high redshift. This may be compared to the actual 
numbers observed, 25 and 14 respectively, and can be seen to differ 
by more than the la of \/A' expected from Poisson statistics. 

The id and 2d marginalised distributions given all three data 
sets are also shown in Figs. |3| and Q It can be seen that adding in 
a third data set tightens the constraints further and the 68 per cent 
confidence limits obtained from these Id distributions are presented 
in Table 1 . 



5 DISCUSSION 

The non-trivial and reassuring result is that the marginalised distri- 
butions of the parameters are consistent when single data sets, or 
any combination of two data sets, are included. This means that no 
single experiment is dominating or pushing the favoured parameter 
values away from the otherwise preferred ones. It is thus tempting 
to suggest that systematic uncertainties in the individual data sets 
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Figure 3. The Id marginalised distributions using (a) pairs of data sets and (b) using all three data sets. The dotted lines denote 68 per cent confidence limits. 



0.8 

0.7 
0.6 
0.5 
0.4 
0.3 




0.2 0.4 0.6 0.1 

Q. 



1.2 r 
1 

coO.8 

D 

0.6 
0.4 




0.2 0.4 0.6 O.i 

Q. 



1.2 
1 

coO.8 
0.6 
0.4 




0.4 0.6 
h 



0.8 



Figure 4. The likelihood function from combining all three data sets as a function of pairs of parameters, in each case marginalised over the remaining two 
parameters. The contours denote the 68, 90, 95 and 99 per cent confidence regions. 



may not be greatly affecting the parameter values that have been 
inferred by this analysis. 

The strong agreement between pairs of data sets can be un- 
derstood qualitatively as follows. The cluster abundance data pro- 
vide a strong constraint on Og; the CMB provides a much weaker 
limit on Og but one which agrees with the cluster abundance limit. 
The IRAS data provides limits on the combination CTg^J/jMS and 
so putting this together with the above constraint on ag produces 
a small allowed range of ^iras- The cluster abundance data also 
favour just a small range in flm. The different degeneracy direc- 
tions in the Q.m,h plane from each of the CMB and IRAS data pro- 
duce a slightly weaker limit on Q.^, but one which agrees with the 
cluster abundance result. Putting together the strong cluster abun- 
dance constraint on flm with either the CMB or IRAS limits in the 
h, iim plane produces a constraint on h; alternatively, just the CMB 
and IRAS data can be used to determine an allowed range in h; and 
each of these routes to h produces the same answer. 

Our results agree well with other recent work which combines 
data sets to set limits on cosmological parameters (Bond & Jaffe 
1998, Gawiser & Silk 1998, Lineweaver 1998, White 1998, Gar- 
navich et al. 1998, Eisenstein et al. 1998 and Roos & Harun-or- 
Rashid 1999). 

The recent Type la supernova results of Perlmutter et al. 
(1998) imply that for a flat universe ilm = 0.28+oog {la statisti- 
cal errors), in good agreement with the values derived in this letter. 
Similar results are also found by Riess et al. (1998). From gravi- 
tational lensing measurements, Falco, Kochanek & Munoz (1998) 



find that Hm > 0.38 at 2 sigma, although Chiba and Yoshi (1999) 
find that a flat universe with = 0.3^q j is preferred. Our value 
for the Hubble constant, h = 0.54, falls at the low end of the latest 
distance ladder measurements (Freedman et al. 1998). The optimal 
baryon density (given flb/i^ = 0.019) is fib = 0.066, and dividing 
by our best fit value of £1^ = 0.36 our optimal baryon fraction is 
0.18 This may be compared to cluster baryon fraction estimates, 
for example. White et al. (1993) find nb/n,n > 0.015 + 0.08/j"^/2. 
Inserting our best fit value for h yields ilb/^m ^ 0.22. Our value 
for the combination OgflU;^ = 0.40 is lower than the one derived 
from measurements from the peculiar velocity field OgiljJ;^ ~ 0.8 
(Freudling et al. 1998). It is closer to the value recently derived 
from cluster peculiar velocities by Watkins (1997) of Cgilj],^ = 



0.44 



h0.19 



Reassuringly, the age of the universe we find (16.7 Gyr) 



is greater than the age of the oldest stars (Chaboyer et al. 1998). 
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